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ABSTRACT 

We discuss the properties of subhalos in cluster-size halos , using a high-resolution statistical sample: 
the Rhapsody simulations introduced in Wu et al. (2012 1. We demonstrate that the criteria applied 
to select subhalos have significant impact on the interred properties of the sample, including the scatter 
in the number of subhalos, the correlation between the subhalo number and formation time, and the 
shape of subhalos' spatial distribution and velocity structure. We find that the number of subhalos, 
when selected using the peak maximum circular velocity in their histories (a property expected to be 
closely related to the galaxy luminosity), is uncorrelated with the formation time of the main halo. 
This is in contrast to the previously reported correlation from studies where subhalos are selected by 
the current maximum circular velocity; we show that this difference is a result of the tidal stripping of 
the subhalos. We also find that the dominance of the main halo and the subhalo mass fraction are 
strongly correlated with halo concentration and formation history. These correlations are important to 
take into account when interpreting results from cluster samples selected with different criteria. Our 
sample also includes a fossil cluster, which is presented separately and placed in the context of the rest 
of the sample. 

Keywords: cosmology: theory — dark matter — galaxies: clusters: general — galaxies: halos - 
methods: N-body simulations 



1. INTRODUCTION 

The abundance and spatial distribution of galaxy clus- 
ters in the universe have played an essential role in de- 
termining cosmological parameters. These properties 
are sensitive to cosmic expansion and the large-scale 
structure growth rate, making clust ers complementar y to 
other cosmological pro bes (see, e.g., Allen et al.|2011 and 
|Weinberg et al.||2012| for recent reviews, and references 
therein). Among multi- wavelength cluster surveys, opti- 
cal surveys provide the largest statistical power in terms 
of the number of identified clusters; this number will dra- 
matically increase with the next generation of wide area 
surveys, including PanSTARRS^ DES^J EuclicjfJ and 
LSSTrj However, the precision cosmology that can poten- 
tially be achieved will be limited by the systematic effects 
invol ved, including clu ster identification and centering 
(e 



relevant observations. A common procedure is to use dark 
matter-only N-body simulations to predict the distribu- 
tion of dark matter particles and halos, and then relate 
halos and subhalos to the o bserved galaxy clusters and 
their memb er galaxies (e.g., Kravtsov et al.|[2004| |Zheng 
et al.|[2005| ). However, resolving subhalos in cluster-size 



halos comparable to the observable limits and associating 
them with galaxies presents additional challenges. As de- 
scribed by the hierarchical structure formation paradigm, 
subhalos accrete onto the main halo through numerous 
merger events and have been substantially influenced by 
the deep gravitational potential of the mai n halo (e.g., 



Ghigna et aLfll998| [ri 

simulating these sub" 



simulating 
resolution (e.g 



Moore et al. 1998 1999). Therefore 



halos require s high mass and force 
Klypin e t al. | 1999[ ) , improved halo finding 
(e.g., Onions et al.||2U12[ ), as well as a ca reful modeling 



Rykoffet al.|20i2t , the nor malization and scatter o f ofHie associated satellite galaxies (e.g., |Reddick et al. 



the richness-mass relation (e.g 
orientation and projection effects" 



Rozo et al. 2009, 2011 
Uohn et al.||2Qi 



(e.g. 



White et al.||2010 Erickson et al.|20TT|), cross-comparison 



with multi-wavelength data ( |Rozo et al.||2012 ), as well as 



uncertainties in theoretical calibrations of halo statistics 
( Wu et al.|[20T0 i. 

To quantify the systematic effects inherent in the mea- 
surements of galaxy clusters, it is essential to generate 
a simulated sample of clusters that is comparable to the 

The Panoramic Survey Telescope & Rapid Response System; 
http: //pan-starrs . if a.hawaii . edu/ 

2 The Dark Energy Survey; http : / /www . darkenergy survey . org/ 

3 http://sci.esa.int/euclid/ 

4 The Large Synoptic Survey Telescope; http://www.lsst.org/ 



To characterize the galaxy populations in galaxy clus- 
ters obtained from deep wide surveys, it is necessary 
to simulate clusters with high resolution (to resolve the 
galaxy content to observable limits) and in a large cos- 
mological volume (to obtain a statistical sample), which 
is computational ly challenging. As discussed in Paper I 
(Wu et al. 2012[ ), in order to achieve large sample size 
and high resolution simultaneously, we have repeatedly 
applied a "zoom-in" or multi-resolution simulation tech- 
nique to develop the Rhapsody sample, which currently 
includes 96 halos of mass lO 14 - 8±O O5 ft,- 1 M , selected from 
a cosmological volume of side length 1 /i -1 Gpc and re- 
simulated with mass resolution 1.3 x lQ s h~ 1 M ( T ) . This 
sample is currently unique in terms of its sample size and 
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resolution and occupies a new statistical regime of cluster 
simulations (see Figure 1 in Paper I). In this second paper, 
we focus on the subhalo population of the Rhapsody 
clusters and give particular attention to the impact of 
formation history and tidal stripping on subhalos. 

The impact of formation history on the observable prop- 
erties of clusters (e.g., galaxy number and distribution) 
is important because it can provide extra information 
or, if not correctly taken into account, introduce bias in 
cluster mass calibration and in cosmological constraints 
that depe nd on the propert ies of galaxies in clusters. For 
example, |Wu et a~ (2008) have shown that if the rich- 
ness of a cluster (the number of galaxies in a cluster 
under a certain selection criterion, used as a cluster-mass 
indicator) is correlated with its formation time, then 
richness-selected clusters will be impacted by assembly 
bias (i.e., early forming halos have a higher halo bias), 
which will in turn impact cluster mass self-calibration 
and cause systematic errors in the inferred cosmological 
parameters. Cosmological studies which use information 
about the halo occupation of galaxies (the number of 
galaxies inside a halo for a given halo mass) also depend 
on an understanding of whether this occupation depends 
on properties other than halo mass. Therefore, it is im- 
perative to characterize these correlations with higher 
precision using a statistical sample relevant for current 
and future surveys. 

In this paper, we discuss how formation history impacts 
the subhalo abundance, subhalo mass fraction, and the 
dominance of the main halo over its subhalos. In par- 
ticular, we focus on the influence of the specific criteria 
used to select subhalos from simulations on the inferred 
properties of the subhalo population thus obtained. One 
of our main findings is that the correlation between sub- 
halo number and formation time sensitively depends on 
this selection criterion. If we use a subhalo selection cri- 
terion that is insensitive to the stripping of dark matter 
particles, the subhalo number and halo formation time 
are not correlated. This result implies that both the halo 
occupation and the cluster richness for halos of a given 
mass are not likely to correlate with formation time. 

This paper is organized as follows. In Sj2] we briefly 
summarize our simulations and halo catalogs, as well as 
the various subhalo selection criteria that we consider. In 

fwe present the statistics of subhalos and the shape 
their spatial distribution and velocity ellipsoid. In <|4j 
we discuss how formation time impacts observational sig- 
natures of subhalos, including the subhalo mass fraction 
and the dominance of the main halo. In fj5] we focus on 
the impact of halo formation time and tidal stripping on 
the number of subhalos. We conclude in [J6j 

2. HALO CATALOGS 

The Rhapsody sample includes 96 cluster-size halos 
of mass M vir = 1O 14 - 8±O O5 /i~ 1 M , re-simulated from a 
cosmological volume of 1 ft~ 3 Gpc 3 . Each halo has been 
simulated at two resolutions: 1.3 x 10 8 /i _1 M Q (equivalent 
to 8192 3 particles in this volume), which we refer to 
as "Rhapsody 8K" or simply "Rhapsody"; and 1.0 x 
1O 9 /i _1 M (equivalent to 4096 3 particles in this volume), 
which we refer to as "Rhapsody 4K." The simulation 
parameters are summarized in Table 1 of Paper I. 

All simulations in this work are based on a ACDM 
cosmology with density parameters fl m = 0.25, Q\ — 



0.75, Qb — 0.04, spectral index n s 
<j g = 0.8, and Hubble parameter h ■■ 



= 1, normalization 
0.7. 



2.1. The Simulations 

The implementation of Rhapsody can be summarized 
as follows: 

1. Selecting the re-simulation targets: We start from 
one of the 1 /i _1 Gpc volumes (named "Carmen") 
from the LasDamas suite of simulation^] and select 
halos in a narrow mass bin io 14 - 8±0 05 ft _1 M0. 

2. Generating initial conditions: We use t he multi 
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scale initia l condition generator Music (Hahn & 
Abel 2011) to generate "zoom" initial conditions 
for each cluster with the second-order Lagrangian 
perturbation theory. 

Performing gravitational evolution: We compute 
the non-linear evolution of each cluster down to z — 



usin g the public version of Gadget-2 (Springel 



2005) 



Identifying halos and subhalos: We u se the adaptive 
phase- space halo finder Rockstar (Behroozi et al. 
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2011a) to assemble catalogs of halos and subhalos 
at 200 output times. Rockstar achieves a particu- 
larly high completeness of the subhalo sample. 

Constructing merger trees: We use th e gravita- 
tional l y consi stent merger tree code by Behroozi 



et al. (2011b) to construct merger trees from the 
halo/ subhalo catalogs. 

We kindly refer the reader to § 2 of Paper I for more details 
on the simulations, the halo and subhalo identification 
and merger tree generation, as well as the mean values 
and variances of the various key properties of the main 
cluster halos in Rhapsody (given in Table 2 of Paper I). 

2.2. Subhalo selection methods 

Subhalos in cluster-size halos are expected host the 
observed satellite galaxies in clusters. For each main halo 
in Rhapsody, we consider the subhalos within its virial 
radius, R v r (based on the spheric al overdensity calcu- 
Bryan fc Norma n||l9 98| A v 



lated with 



i94c 



at z=0). 



We characterize each subhalo by the maximum circular 
velocity of the dark matter particles associated with it, 
Vmax, defined at the radius r = r max that maximizes 
y/GM(< rj/r: 



l GM{< r n 



(1) 



This quantity is often used as a proxy for subhalo mass, 
since the mass of a subhalo itself is typically not a well 
defined quantity in the simulations due to ambiguities 
about how to separate subhalos from the background 
density of the main halo. We focus on v max at two 
different epochs during the evolution history of a subhalo: 

• v : the value of v max measured at z = 0, a quantity 
related to the current subhalo mass. 

5 http : //lss .phy . vanderbilt . edu/lasdamas/ 
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• v p k: the highest v max value in a subhalo's history, 
a quantity related to the highest subhalo mass in 
its entire history. 

Since subhalos experience strong tidal stripping after their 
accretion onto the main halo, the two are not identical. 
The parameter v p k is more closely related to the luminos- 
ity and stellar mass of satellite galaxies than vo, because 
the stellar component of a galaxy is denser and less easily 
stripped than the more extended dark matter component. 
Even though a halo could lose dark matter particles at 
its outskirts, the galaxy in its core can remain intact 
for a longer time. In f act, the satellite gal axy may even 
continue to grow (e.g., Wetzel et al.||2012 ). Therefore, a 
quantity that is unaffected by stripping is expected to 
provide a better proxy for the stellar mass of galaxies. 

In particular, subhalo abundance matching models 
based on properties that are less impacted by stripping 
have bee n shown to better ag r ee wit h obs ervations. For 
examp le, Nagai & Kravtsov (2005) and Conroy et al. 
(20061 have shown that using v ac (i.e., v max at the time of 
accretion of the subhalo) when selecting subhal os better 
reproduces t he statistics of observed galaxies. |Reddick| 
et al. ( 20121) have further demonstrated that an abun- 
dance matching model based on v p k provides a better 
fit than either vo or v ac to the galaxy two-point correla- 
tion function and the conditional stellar mass function 
for galaxies in groups from the Sloan Digital Sky Survey 
(SDSS). 

In addition to v max , we also investigate the mass of 
subhalos at two different epochs: 

• Mq\ the current subhalo mass, defined by the par- 
ticles bound to the subhalo according to the imple- 
mentation of ROCKSTAR. 

• Afpk: the highest mass in a subhalo's entire assem- 
bly history. 

In general, the difference between Mq and M p k is analo- 
gous to the difference between vo and v p k, except that the 
circular velocity is less affected by stripping than mass is. 

It is important to note here that finite resolution in N- 
body simulations l eads to the "overmerging" effect (e.g., 
Klypin et al.|1999 |: small subhalos tend to fall below the 
resolution limit before they merge with the central object. 
However, the stellar component associated with subhalos 
is expected to survive longer than the simulated subhalos. 
A detailed discussion of the resolution dependence of 
this effect and the associated completeness limits of the 
subhalo populations will be presented in a separate paper 
(Wu et al., in preparation). 

3. SUBHALO STATISTICS AND DISTRIBUTIONS AT Z = 

In this section, we focus on the statistical properties 
of the subhalos in Rhapsody, as well as the shape of 
the subhalos' spatial distribution and velocity ellipsoid. 
In particular, we explore the impact of resolution and of 
the various selection criteria described in the previous 
section. 

3.1. Subhalo mass function 

The mass function of subhalos has been shown to follow 
a power law for low mass subh alos and an exponential 
cutoff for massive subhalos (e.g., Gao et aL]|2004| |Angulo| 



et al. 2009 Giocoli et al. 2010). In this section, we 
investigate the validity of this form when using the various 
criteria for subhalo selection . 

Figure [l] shows the number of subhalos above a given 
threshold of v (left), v p k (middle), and M (right). The 
blue/red curves correspond to subhalos within i? v ; r of 
RHAPSODY 8K/4K halos, while the transparent grey 
curves correspond to individual halos in Rhapsody 8K. 
The black curves correspond to all halos and subhalos 
within 7 /i _1 Mpc around the center of the main halo, a 
region where re-simulated halos are well resolvecj^J The 
blue dashed lines indicate the best-fit power-law slopes of 
the distribution functions of the 8K sample. We find that 
the slope of our Mq function is slightly shallowe r than|De 



Lucia et al. (2004), Boylan-Kolchin et al. (20101, and |liao 
et al.| ( |20P2 ), which is plausibly attributed to the different 
mass definition. At the same time, our v function is in 



good agree ment with the resu lts of |Boylan-Kolchin et al.| 
d2010| ) and |Wang et al.| ( |2012h based on Milky Way-size 
laios! The large halo-to-halo scatter shown here has 
implications for comparisons of satellite statistics in the 
Milky Way with simulations, given that thus far these 
comparisons have been done with a small number of simu- 
lated halos. Scatter in the properties of satellites between 
halos likely reduces the current tension with observations 
of massive dwarf galaxies (e.g., Purcell & Zentner 2012). 
A statistical sample for galactic subhalos from simula- 
tions, as well as a larger sample of observed systems, is 
required to verify these results. 

The 4K and 8K subhalo mass functions deviate from 
a power law at different values of vo and v p k, clearly 
indicating the dependence of the completeness limit for all 
mass proxies on resolution. The gray curves demonstrate 
the significant scatter in subhalo abundance from halo 
to halo in our sample. Given this large scatter, we next 
explore how well the number of subhalos in different bins 
are correlated, i.e., how sensitive a richness estimator 
would be to a different selection threshold. In Figure [2] 
we assign subhalos into bins of 50 km s -1 using vo (left) or 
v p k (right), starting from 50 km s _1 (we note that the first 
bin in either case is incomplete, and the last bin includes 
all subhalos beyond 250 km s _1 ). We compare each pair 
of bins and find that the subhalo counts are only weakly or 
moderately correlated between bin^H This indicates that 
a halo that is rich in massive subhalos is not necessarily 
also rich in low-mass subhalos. If the satellite galaxy 
populations in clusters in different luminosity bins follow 
the statistical distribution shown here and have such low 
covariances, they could potentially provide independent 
information for mass calibration. 

3.2. Scatter of subhalo number 



Boylan-Kolchin et al. (2010) and Busha et al. 
have shown that for galactic halos (M < 10 1,s b h 



(2011) 



Mr. 



6 To decide the size of the well-resolved ambient region around the 
main halo, we compare halos in the re-simulated region (composed 
of high-resolution particles only) with those in the original Carmen 
simulation. At 7 h~ Mpc, the re-simulated regions recover the halo 
population in the corresponding region in the Carmen simulation. 
In addition, at 7 /i Mpc, the number of low-resolution particles is 
less than 4%. although it varies from halo to halo and is sometimes 
0%. 

7 Throughout this work, we use rank correlation, which makes 
our results insensitive to outliers. 
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Figure 1. The number of subhalos above a given threshold vo (left), v p k (middle), and Mq (right). The blue/red curves represent the 
mean number of subhalos within the i? v ir of RHAPSODY 8K/4K halos, while the black curves represent all halos and subhalos within 7 
h~ Mpc around the center of the main halo in the re-simulation. The thin gray curves in the background show the subhalos for individual 
Rhapsody 8K halos. The blue dashed line indicates the slope of the 8K sample in the regime where the subhalo sample is complete. 
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Figure 2. Correlation between the number of subhalos in different bins. Subhalos are binned by vo (left) or Vpk (right), with bin size 50 
km s — 1 , starting from 50 km s" 1 . The subhalo numbers between different bins have moderate or weak correlation. 



the distribution function of the number of subhalos (N) 
deviates from the Poisson distribution when (N) is large. 
We repeat this analysis for our sample of halos of sig- 
nificantly higher mass. Figure [3] shows the scatter of N 
under different subhalo selection methods and thresholds. 
The left panel shows the ratio between the measured scat- 
ter in the sample a = ^/Var[N] and the Poisson scatter 

cpoisson = VT^O- The right panel presents the second 
moment of the subhalo number distribution 



^/(N(N-l)) 
(TV) 



(a = 1 for Poisson) 



(2) 



Both quantities are measures of how the distribution 
deviates from the Poisson distribution. 

In both panels, the a;-axis corresponds to (N) for a given 
selection threshold, allowing direct comparison between 



different subhalo selections. Each curve corresponds to 
a different selection method, and several corresponding 
thresholds are marked on each curve (in units of km s _1 
for maximum circular velocity and log 10 ^"'Mg for mass). 
For low thresholds or high (N), the scatter deviates sig- 
nificantly from the Poisson scatter. In addition, a has 
a trend similar to that of g/gp j s son and only slightly 
deviat es from unity. As discussed in |Boylan-Kolchin et al.| 
(2010), slight deviations of a from unity can correspond 
to large deviations from the Poisson distribution. We also 
note that the additive boost of the variance has a similar 
trend as a, since (a 2 - (N))/(N) 2 = a 2 - 1. 

From Figure [3j our results for different selection criteria 
can be summarized as follows: 

• vo vs. Vpk (red vs. blue): v p k selection gives less 
scatter and is closer to the Poisson distribution, 
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Figure 3. Scatter of the number of subhalos for various selection methods. The left panel corresponds to the ratio between the sample 
scatter and the Poisson scatter, while the right panel corresponds to the second moment of the distribution. The x-axis is the mean number 
of subhalos for each selection method, and several of the thresholds are marked on each curve. Here we compare various cases: (1) vo 
vs. Vpk (red vs. blue), (2) Mo vs. Mpk (purple vs. green), (3) 4k vs. 8K (purple dashed vs. purple solid). In each pair, the former selection 
method leads to extra non-Poisson scatter. This trend indicates that both stripping of subhalos and insufficient resolution can induce extra 
non-Poisson scatter and might lead to a scatter greater than the values in observations. For v p j< selection, a = 1.005 for sufficiently large 
(N). 

indicating that stripping of subhalos can introduce The shape parameters are defined through the distri- 

cxtra non-Poisson scatter. bution tensor: 



• M vs. M pk (purple vs. green): M pk selection is 
closer to the Poisson distribution, which can also 
be understood with stripping. 

• v p k vs. Mpk (blue vs. green): similar. Both prop- 
erties are computed before a subhalo's infall and 
behave similarly. 

• Mq 4K vs. 8K (purple dashed vs. purple solid): 
8K is closer to the Poisson distribution, indicating 
that insufficient resolution can introduce extra non- 
Poisson scatter. 

In all cases presented here, the v p k selection provides 
the smallest non-Poisson scatter, with an asymptotic 
value of a = 1.005 for sufficiently large (N). Our results 
show that stripping and insufficient resolution can lead 
to extra non-Poisson scatter. This trend can also explain 
the modest difference between our values of non-Poisson 



( (n -( ri )) (rj .-( r .))) 



(3) 



scatter and those quoted in Boylan-Kolchin et al. (2010) 
who have higher resolution but used subhalo mass. 



3.3. Subhalo spatial distribution and kinematics 

In Paper I, we have discussed the shape and velocity 
ellipsoid of dark matter particles of the main halo. It 
is interesting to see how closely the subhalos follow the 
distribution of dark matter in position and velocity space, 
where differences exist, and how these depend on the 
specific selection of subhalos. We thus present analogous 
measurements for subhalos, which are selected with vo 
and v pk . 



where n is the i th component of the position vector r of 
a subhalo. The eigenvalues of ly are sorted as Ai > A2 > 
A3, and the shape parameters are defined as: a = y/Xx, 
b = \/A2, c = \A3- We present the dimensionless ratios 
b/a and c/a. In addition, the triaxiality parameter is 
defined as 

a 2 -b 2 

T = -2 2 ■ 4 

T rj I (a > b ft c) indicates a prolate halo, while T f=a 
(a ft b > c) indicates an oblate halo. Intermediate values 
of T correspond to triaxial halos. 

Analogously, th e velocity ellipsoid is defined as (e.g., 
White et al.|[2010] ): 



((Vi - (Vi))(Vj - (Vj))) 



(5) 



where is the i component of the velocity vector. 
Sorting the eigenvalues of the velocity ellipsoid as Ai > 
A 2 > A3, one can again define = \fX\, b^ v ' = VA2, 

The scatter of the velocity dispersion along different 
lines of sight can be calculated from the eigenvalues of 
the velocity ellipsoid tensor as 

(^L) = ^(Ai + A 2 + A 3 ) (6) 

(^L) 2 = ^(a? + a2 + a2-a 1 a 2 -a 2 a 3 -a 3 a 1 ) . (7) 
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Figure 4. The distribution of the shape and velocity ellipsoid parameters for dark matter particles and subhalos within R V1I of the main 
halos. The upper panels correspond to the shapes of spatial distribution, and the lower panels correspond to the shapes of velocity ellipsoid. 
As shown in Panels (a) and (b), dark matter particles tend to be more prolate than subhalos. Panel (c) shows that the velocity ellipsoids of 
subhalos tend to be more elongated. Panel (d) shows that the statistical error of line-of-sight velocity dispersion measurements is larger for 
subhalos than for dark matter particles. 
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We use subhalos within i? v j r selected above a given thresh- 
old without imposing weighting. 

Figure [4] shows the cumulative distribution function 
of the shape and velocity ellipsoid parameters of dark 
matter particles (solid), as well as of subhalos selected 
with v > 100 km s _1 (dotted) and subhalos selected 
with v p k > 150 km s _1 (dashed). The upper panels 
show the axis ratios and triaxiality parameters. Subhalos 
selected by vo have a distribution closer to spherical than 
both Vpk-selected subhalos and dark matter particles. In 
addition, the dark matter distribution tends to be more 
prolate than the subhalo distribution. 

The lower panels of FigureHJshow the shape parameters 
of the velocity ellipsoid ana the line-of-sight scatter of 
the velocity dispersions. We find that the velocities of 
subhalos tend to be more anisotropic than those of the 
dark matter particles (based on c^/a^"'). In addition, 
Saf os is higher for subhalos than for dark matter particles. 
This trend can be understood based on the findings of 
| White et al.| ( |2010[ ) : the motions of subhalos tend to be 
anisotropic because they can retain their infall velocities 
for a long time. In contrast, the merged and stripped 
material is dynamically older and contributes to a well- 
mixed isotropic velocity distribution. 

We summarize several trends in Figure [4] and propose 
explanations: 

• The distribution of subhalos tends to be more spher- 
ical than that of dark matter. This can be under- 
stood by the fact that subhalos have a shorter re- 
laxation time than dark matter particles, because 
treiax ~ (R / v) (N / la N) and the N for subhalos is 
much smaller. 

• The velocity ellipsoid of subhalos tends to be more 
elliptical than that of dark matter. This can be 
explained by the anisotropic motions of subhalos 
accreted as a group (as discussed above). 

• Subhalos selected with vo tend to have a more 
spherical distribution than those selected with Vp^. 

Similar trend exists for /a^ but does not ex- 
ist for b^/a {v K This trend can be understood 
through stripping: a v p k selection tends to include 
more highly-stripped subhalos than v , and highly- 
stripped subhalos tend to be on more elliptical or- 
bits (stripping will be stronger for those orbits that 
have a smaller pericenter), leading to the higher 
ellipticity measured with v p k. In addition, we note 
that overmerging will also make the distribution 
of subhalos more spherical, because overmerging 
tends to eliminate highly-stripped subhalos, which 
tend to have more elliptical orbits (again due to the 
smaller pericentric distance). 

We note that for both vo and v p k selection, when we 
increase the threshold, the ellipticity slightly increases. 
This trend can be explained by the statistical biases aris- 
ing when a smaller number of subhalos is used to measure 
the ellipsoids. To confirm this, we randomly select a num- 
ber of subhalos within i? v j r with (N) matching either the 
number of subhalos obtained with the vo- or v p k-selection, 
and we recover the trend that a smaller number of subha- 
los always leads to a higher inferred ellipticity. Although 



it is possible that large radius is weighted more in the 
position tensor and that small radius is weighted more in 
the velocity tensor, we find that the radial distribution of 
massive subhalos alone is insufficient to explain the trend 
with selection threshold. 

In addition, we observe that the difference between 
shapes measured by dark matter particles and subhalos 
has a slight trend with the formation history — for halos 
that experienced recent major mergers, the subhalo dis- 
tribution tends to be much rounder than the total dark 
matter particle distribution. This could also be related 
to the fact that subhalos have a shorter relaxation time 
than dark matter particles. 

Finally, we note that in all cases, the differences in 
axis ratios are at the level of a few percent. Observa- 
tionally, these differences are likely to be overwhelmed 
by the scatter due to line-of-sight projection, viewing 
angle, spectroscopic sample selection, etc. These effects 
are likely to depend on the environment and details of the 
observational techniques employed. It would be interest- 
ing to investigate whether and, if so, how the difference 
between the velocity ellipsoids for the var ious selection 
criteria depends on environment (see also Faltenbacher 
|2010[ who have demonstrated a dependence ot subhalo 
kinematics on environment in the Millennium Simula- 
tion). An analysis of the environmental dependence of 
halo properties will be deferred to a future paper. 

4. CORRELATION BETWEEN SUBHALO PROPERTIES AND 
FORMATION HISTORY 

In Paper I, we have focused on the impact of the for- 
mation history on the density profile and the phase-space 
structure of the cluster halos. Here, we investigate the im- 
pact of the formation history on the subhalo population. 
In Figure [5l we present the correlation of eight quantities 
measured from our sample: four for subhalo properties, 
three for formation time, and one for the halo profile. 

4.1. Formation history parameters 

As discussed in Paper I, to quantify the formation 
history of halos, we adopt an exponential-plus-pow er law 
model with two parameters ( McBride et al.|[2009 ) 



M(z)=M (l - 
d\nM 

— »7"P 



zfe 



when z << 1 



(8) 
(9) 



Thus, 7-/3 provides a measure of the late-time accretion 
rate. 

We use Z1/2, the earliest redshift that a halo obtains 
half of its mass, as the formation time proxy throughout 
the paper, but note that using formation time proxies 
based on fitting functions (as those studied in Paper I) 
lead to similar results. In addition, we include the redshift 
of the last major merger of each main halo, Zi mm . 

We note that the formation history is directly reflected 
in the halo density profile, which has been studied in 
detail in Paper I. For the completeness of our correlation 
analysis, we also include the halo concentratio n defined 
for the Na varro-Frenk-White (NFW) profile (Navarro 
eTaLlp97| ): 



P{r) = 5 C 

Pent (r/r s )(l + r/r s )' 



(10) 
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for which the concentration parameter is defined as 

Rvir 

cnfw = • 



(11) 



We kindly refer the reader to Paper I for details on the 
fitting procedure. 

4.2. Subhalo mass fraction 

We next investigate the mass of the main halos that is 
contained in subhalos. The subhalo mass fraction can be 
defined as 



fsub(Mth) 



1 



E 



(12) 



M sub >M th 



and can be used as an indication of a recent major merger 
event. For example, if a massive subhalo accreted onto 
the main halo only recently, it retains most of its mass 
and contributes to a large f su b- We show below that 
f su b is correlated with halo formation time and can thus 
be used as an indicator of the state of relaxedness o f 
the halo (e.g., |De Lucia et al ||2004| |Shaw et al.||2006 ) 



Here we choose M t h = 10 w Ii~ 1 Mq, which approximately 
corresponds to our completeness limit (see Fig. [I]), but 
note that the correlations presented below are insensitive 
to the specific choice of M t h- 

The third row and column of Figure[5]correspond to f su b 
and show, as expected, that f su b is strongly correlated 
with Z1/2, 7-/3, and Cnfw- That is, halos of higher f su b 
tend to be late-forming, with high late-time accretion 
rates and associated low concentration. 

The subhalo mass fraction itself is also a quantity of 
observational interest. For example, it can be inferred 



from gravitational lensing (e.g., 


Dalai & Kochanek||2002 


Vegetti et al. 2012 


Fadely & Keeton|2U12|). Accurate mod- 


eling ot the subhal 


o mass traction is essential for the study 



of the lensing flux ratios (e.g., Xu et al.]|2009[ ). For the 
relatively massive systems considered here, we find that 
f su b is strongly correlated with the mass of the most mas- 
sive subhalo, despite the fact that this halo contributes on 
average only ~ 20% of the total subhalo fraction. We also 
find that the subhalo fraction is strongly correlated with 
the dominance of the main halo (related to the luminosity 
gap between brightest and second brightest galaxies), as 
well as with formation time and concentration. Because 
strong lensing clusters tend to have higher than average 
halo concentrations, it is important to take this correla- 
tion into account when interpreting measurements of the 
subhalo mass fraction from strong lensing. 

4.3. Mass contributed by merged subhalos 

In the previous subsection, we addressed the mass con- 
tributed by the present subhalo population. We now 
investigate the mass that was brought into the main halo 
by all merging events in a halo's history. Figure [6] shows 
the contribution to the main halo mass from merged 
subha l os. T his has also been explored by, e.g., |Berrier| 
et al. ( |2009[ ) for lower mass systems. The x-axis corre- 
sponds to the ratio of subhalo mass to main halo mass, 
[i = M sub /M™ am , and the y-axis corresponds to the frac- 
tion of main halo mass contributed by subhalos above a 
given fi. 

Here we consider two types of subhalos. The first type 
is those subhalos that have merged into the main halo 



and can no longer be identified; for this type of subhalo, 
we use its mass when it accreted onto the main halo, M ac . 
These subhalos are represented by the blue curve, and the 
blue shaded region corresponds to its standard deviation. 
Our results indicate that, on average, 60% of the main 
halo's mass comes from merged subhalos with [i > 10~ 4 ; 
however, this number varies greatly from halo to halo. 

The second type is those subhalos that still survive 
today (i.e., that can be identified by the halo finder at 
z = 0). For this type, we also use M ac . Since subhalos 
tend to lose a significant amount of mass due to tidal 
stripping inside the main halos, using M ac ensures that 
we include the mass that once belonged to subhalos but 
later got stripped by the main halo. We do not explicitly 
count the subhalos that merge into other subhalos because 
their masses have already been included in the surviving 
subhalos. The red curve corresponds to the sum of the 
first and the second types. We find that more than 90% 
of the mass can be attributed to halos with /i > 10~ 4 
that were accreted onto the main halo. 

Finally, the black curve shows the contribution to main 
halo from subhalos that survive at z = 0. We use the cur- 
rent mass of these subhalos, M ; this quantity is equiva- 
lent to the subhalo mass fraction for different mass thresh- 
olds, fsub(> m)- Subhalos with p, > 10~ 4 constitute a 
little less than 20% of the mass of the main halo. 

4.4. Dominance of the main halo 

In this section, we study the difference between each 
main halo and its largest subhalo, a quantity motivated 
by the definition of the so-called "fossil groups" . Obser- 
vationally, these systems are defined as having a large 
magnitude gap between the brightest and second bright- 
est galaxies, in addition to being X-ray luminous (e.g., 
Tremaine fc Richstone||T977| | Jones et al?] [20031 |Miller| 
et al. 2012 1 Fossil systems are often interpreted as a 
population of galaxy groups that have assembled at early 
times and have not undergone a recent major merger. 

Since predicting the optical and X-ray properties of 
our halos is beyond the scope of this paper, we define a 
related property, the ratio of v p k for the main and the 
first subhalo (the subhalo with highest v p k): 



D = 



"pie 



-1st sub ' 



(13) 



We note that using the second subhalo leads to the same 
trend presented below. 

In Figure |5j the fourth column and row show that D is 
correlated with the formation time, late-time accretion, 
concentration, and subhalo mass fraction. We also note 
that all our main halos have similar vp^f 111 ; therefore, the 

scatter in D is almost completely determined by v^* sub . 

The trends observed in Figure [5] can be understood 
as follows: Since Vpf.* sub indicates the maximum of the 
subhalo mass that accretes onto a main halo, a main halo 
with a low v 4 ^ sub has fewer massive subhalos accreting 
onto it (this is also reflected by its low f su b). With 
relatively fewer incoming subhalos, to achieve the same 
mass today, these halos must have obtained most of their 
mass at early times and have undergone slow accretion 
at late time, thus leading to the low 7-/3 and the high 
concentration. 
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Figure 5. Correlation between the subhalo properties and the formation time proxies. The top two rows show the number of subhalos 
selected with vo and v p i;. In each case, the threshold is chosen to have an average of 100 subhalos per main halo. Correlations tend to be 
weaker with subhalos selected by v p ^ than by vq, except for z; mm . 
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Figure 6. The contribution of mass to the main halo from subhalos 
above a certain mass ratio ft = M aub /M™'" n . The blue curve 
corresponds to the contribution from subhalos that have merged 
with the main halo, and we use the subhalo mass at accretion, 
M ac ; the shaded area corresponds to the standard deviation of the 
sample. The red curve includes merged subhalos and those that are 
surviving, for which we also use M ac . The black curve corresponds 
to subhalos that are still surviving today, for which Mq is used, and 
is equivalent to f 3U b(> AO- 
While we were prepari ng this manuscript, we learned 
about the related work of Hearin et al. (2012), who have 
studied the "magnitude gap" of the two brightest cluster 
members, which is analogous to our dominance parameter 
D. These authors have found that for SDSS groups of a 
given velocity dispersion, clusters with high magnitude 
gap tend to have low richness, and this correlation can 
in turn reduce the scatter in the mass inferences using 
optical mass tracers. In contrast, we find that the number 
of subhalos selected with v p k is not correlated wit h D for 
halos of the same mass. Since the results from IHearinl 



et al. (2012) are based on velocity dispersion rather than 
mass, a fair comparison between our results and theirs 
will require further consideration of the scatter in velocity 
dispersion, scatter of galaxy luminosity at a given v p k, 
as well as observational selections, which are beyond the 
scope of the current work. 

4.5. The curious case of Halo 572: An outlier and a 
fossil cluster 

In the Rhapsody sample, we find one peculiar halo — 
Halo 572 — which is a prominent outlier in formation 
time (highest ^1/2) and occupies the tail of many halo 
properties as well as the corner of several scatter plots 
(marked as stars in Figure [5]). It has unusually high cnfw 
(2.7a deviation from the mean) and central dominance 
(3.2cr). It also has one of the lowest late-time accretion 
rates 7-/3 (2.7cr), / su & (1.6er), and subhalo numbers 
selected with several different criteria. This halo obtained 
most of its mass at early time and nearly stopped accreting 
mass at late time, leading to these extreme properties. 
Images of the evolution of Halo 572 are shown in Figure 3 



of Paper I, where it is evident that this halo had an 
atypical formation history. We find that Halo 572 does 
not live in an atypical environment on large scales. 

The high central dominance indicates that, if such a 
halo is observed, it would likely have a large luminosity 
gap between the brightest and second brightest galaxies, 
and its high concentration will make it X-ray luminous. 
Therefore, we expect that this halo will host a cluster 
that satisfies the criteria of a "fossil." In addition, it is a 
"real" fossil cluster in the sense that it has an unusually 
early formation history. Studies of fossil groups in both 
simulations and observations have come to a range of 
conclusions, with debate about whether fossil groups have 
distinct assembly histories or are merely an intermediate 



state in gala xy formation (see, e.g., Cui et aLpOll Sales 
et al.|[2007| and references therein). Halo 572 presents 



a case of a distinctively early formation history and the 
consequential properties. From these results we conclude 
that it is highly probable that "real fossils" exist in the 
universe but that they are very rare; thus, they require 
more stringent selection criteria to be distinguished from 
some transient states of cluster formation. 

5. CORRELATION BETWEEN FORMATION TIME AND 
SUBHALO NUMBER: THE IMPACT OF SUBHALO 
SELECTION 

In this section, we investigate in detail the correlation 
between formation tim e and subhalo abundance (as previ- 
ously explored by e.g., [Zentner et al.|2 005 Wcchsle r et al.| 
20061 IGiocoli et ai.||2010|). We show that this correlation 



is mainly caused by subhalo stripping and insufficient 
resolution. However, in the real universe, stripping of 
dark matter particles is less relevant to the galaxy content 
in clusters, and this correlation is therefore not expected 
to exist for observable galaxies. 



Zentner et al. ( 20051) found that early-forming halos 
tend to have fewer subhalos. In their study, subhalos are 
selected with a threshold on vo (v max at z = 0). This 
trend has been explained by the fact that in early-forming 
halos, subhalos tend to accrete at early time and are more 
likely to be destroyed, which leads to a low number of 
subhalos. However, we find that this correlation strongly 
depends on how subhalos are selected. In what follows, 
we explore the dependence of this correlation on various 
subhalo selection criteria and understand the trend by 
investigating the accretion and stripping of subhalos. 

Figure [7] demonstrates how subhalo selection based 
on v or Vpk can lead to different correlations between 
formation time and subhalo number. The left panel 
corresponds to selecting subhalos with vo > 100 km s _1 , 
and the right panel corresponds to v p k > 150 km s _1 . 
We note these two thresholds correspond to roughly the 
same number of subhalos. Each point corresponds to 
a main halo in our sample. The s-axes correspond to 
the number of subhalos under either selection criterion, 
and the j/-axes correspond to the formation time proxy 
z i/2- When subhalos are selected based on vo (left panel), 
formation time and subhalo number are significantly anti- 
correlated. However, when subhalos are selected based 
on Vpk (right panel) , this anti-correlation no longer exists 
in our sample. 

This lack of anti-correlation can be understood as fol- 
lows. When subhalos are selected with a given threshold of 
vo, the stripping of subhalos directly impacts the subhalo 
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Figure 7. The impact of subhalo selection on the correlation between subhalo number and formation time Zi/2- The left panel corresponds 
to selecting subhalos using vo > 100 km s — 1 , while the right panel corresponds to v p k > 150 km s — 1 . Although the vo selection presents 
significant anti-correlation between subhalo number and formation time, the Vpk selection presents no such correlation. This trend can be 
explained by the stripping of subhalos, as demonstrated by the following two figures. 




Figure 8. The evolution of subhalo number, split by quartiles of z x j 2 - Left: subhalos with vo > 100 km s 1 ; right: subhalos with 
v pk > 150km s _1 . When subhalos are selected with vo, the subhalo number of early-forming (blue) and late-forming (red) halos split at 
2 = 0; however, when subhalos are selected with Vpk, there is no clear split of halo number at z = 0. This trend is reflected by the difference 
in the correlation seen in Figure [7] 
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Figure 9. Cumulative distribution function of <5 Vmax = (v p k — 
v o)/ v pk> an indication of the amount of stripping experienced by 
subhalos. We split the main halo by the formation time 21/2 • 
Early-forming halos (blue) tend to have subhalos with higher 5 Vmax 
(stronger stripping) on average than late-forming ones (red). If 
subhalos are selected with vo, highly-stripped subhalos tend to fall 
below the threshold, leading to a low subhalo number. This can 
explain the correlation seen in the vo selection in Figures [7] and [8] 



0.2 1 , , , — 1 — [ 




<N(> V cut or > M cut )> 



Figure 10. Correlation between subhalo number and 21/21 f° r 
different subhalo selection methods. When subhalos are selected 
with vo or Mo (red and purple), an anti-correlation exists for all 
thresholds; when subhalos are selected with v p ^ or Mp^ (blue and 
green), the anti-correlation is greatly reduced or non-existent. In 
addition, the comparison between 8K and 4K (solid and dashed 
of the same color) sample shows that the anti-correlation can be 
enhanced by insufficient resolution. 



number: for a halo that assembled earlier, its subhalos 
experience stripping for a longer time and with a higher 
intensity (because of the high halo concentration), and 
the subhalos' masses and vo tend to be greatly reduced. 
Therefore, fewer subhalos remain above the v threshold. 
On the other hand, when we select subhalos using v p k, the 
stripping of subhalos does not directly impact the subhalo 
number, as long as those subhalos are still identifiable. 

To support our argument above, we investigate the 
evolution of the subhalo population for these two subhalo 
selection criteria. Figure [8] shows subhalo number as a 
function of the scale factor a (we plot in log a to emphasize 
the late-time behavior). As in the previous figure, the 
left /right panel corresponds to the v /v p k selection. In 
both panels, we plot the subhalo number evolution for 
the highest Z1/2 quartile (blue) and lowest Z\i<i quartile 
(red). 

In the left panel of Figure [HI we can see that for early- 
forming halos (blue), the subhalo accretion rate is high 
at early time but suddenly declines after a s» 0.6. On 
the other hand, for the late- forming halos (red), their 
subhalo accretion rate is high at late time. At a ~ 1, the 
early-forming halos have fewer subhalos than late-forming 
halos. In the right panel, although the early-forming halos 
have declined in subhalo accretion rate at late time, their 
subhalo number still grows at late times, and their subhalo 
numbers are similar to late- forming halos at a ~ 1. 

The different trends in both panels can be attributed 
primarily to the stripping of subhalos. Figure [9] shows 
the cumulative distribution of the fractional change of 
v max of subhalos 

Sv _ = V -J^Z^ . (14) 

Vpk 

This quantity can be used as a measure of the amount of 
stripping experienced by subhalos. Higher S VmBX indicates 
that a subhalo has experienced stronger or longer strip- 
ping and has lost more mass. For early-forming halos 
(blue), subhalos on average have higher <5 Vmax , indicating 
that these subhalos experience more stripping and their 
v max is reduced more. As a result, if we select subhalos 
using vo, we tend to exclude subhalos that have experi- 
enced more stripping. These subhalos will however be 
included if we select subhalos using v . 

Therefore, the correlation between formation time and 
subhalo number seen in a selection on vo can be attributed 
to the exclusion of highly-stripped subhalos. Since subha- 
los selected with vo have less observational relevance than 
those selected with v p k, our results imply that cluster 
richness is unlikely to be correlated with the formation 
time of the halo in observations. 

So far we have been using two specific selection thresh- 
olds for v and v p k- Here, we investigate how our results 
depend on the selection threshold. In Figure [TUJ we 
present the correlation between subhalo number and halo 
formation redshift, Cor(N, 2:1/2), where N is the subhalo 
number above some selection thresholds. We discuss four 
selection methods: vo, v p k, Mq, and M p k- For each selec- 
tion threshold, we compute the mean number of subhalos, 
(N) . Using (N) as the 2-axis allows us to put these curves 
on the same figure. 

In Figure [To} the different magnitudes of correlation can 
easily be seen. Here we compare four pairs of selection 
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methods: 

• vo vs. v p k (red vs. blue): the former has stronger 
correlation with Z\/2 due to subhalo stripping, as 
discussed above. 

• M vs. M p k (purple vs. green): the former has 
stronger correlation, for the same reason as above. 

• v vs. v 4K (red solid vs. red dashed): the latter 
has stronger correlation, indicating that an unphys- 
ical correlation can be introduced by insufficient 
resolution. 

• Vpk vs. Vpk 4K (blue solid vs. blue dashed): the 
latter has stronger correlation, indicating that using 
Vpk does not mitigate the impact of resolution. For 
other quantities, comparisons between 8K and 4K 
show the same trend. 

Since v p k is more relevant for observations than vo, Mq, 
and Mpk, the lack of correlation when selecting by v p k 
indicates that, observationally, the formation time of a 
galaxy cluster is unlikely to be inferred from the number 
of galaxies alone. Thus, richness-selected galaxy clusters 
are unlikely to be biased in terms of their formation time, 
implying that the effect of assembly bia s may be negligi ble 
for cluster cosmology self-calibration ( |Wu et _aL 2008). 

We now return to the discussion of Figure [5j We have 
shown that various observables, including the subhalo 
fraction, the central dominance, and the concentration, 
are highly correlated with formation time and the amount 
of late-time accretion. These observables are also corre- 
lated with the number of subhalos selected by the current 
maximum circular velocity, Vo- However, these correla- 
tions largely disappear when selecting subhalos based 
on v p k, which is expected to be more strongly correlated 
with galaxy luminosit y or stellar mass ( Conroy et al.|2006 
Reddick et al. 12012). This reduces the likelihood that 



these observables provide additional mass information for 
richness-selected samples of galaxy clusters. 

6. SUMMARY AND DISCUSSION 

In this paper, we have presented the key properties 
of the subhalo populations in the Rhapsody cluster 
re-simulation project, a sample of 96 halos of M v - lr = 
10 i4.8±o.05^-i Mq; resolved with approximately 5 x 10 6 
particles inside the virial radius. We focus on the effect 
of formation history on the subhalo population. We find 
that this effect depends on subhalo selection criteria and 
resolution, which need to be carefully taken into account 
to make observationally relevant inferences. Our findings 
can be summarized as follows: 

1. Subhalo statistics: In §3.11 we show the subhalo 
mass function for several subhalo mass proxies: v , 
v p k, and M . We find that for a given halo, the 
numbers of large and small subhalos a re o nly mod- 
erately correlated with each other. In §3.2[ we com- 
pare the scatter in subhalo number under different 
selection criteria and resolutions, finding that sub- 
halo stripping and insufficient resolution can lead 
to extra non-Poisson scatter. The least stripped 
proxy, Vpk (8K), still has a small amount of residual 
scatter above Poisson statistics, corresponding to a 
constant value of a = 1.005. 



2. Shape of s patia l distribution and velocity el- 
lipsoid: In £3.3 we compare these quantities mea- 
sured from subhalos selected with vo and v p k, as 
well as from dark matter particles. We find that 
dark matter particles tend to have a more prolate 
distribution than subhalos, and that subhalos show 
a higher line-of-sight scatter of velocity dispersion. 
Subhalos selected with v p k are slightly more ellipti- 
cally distributed than those selected with vq. 

3. Formation history and subhalo properties: 

We have quantified the correlations between various 
subhalo properties and halo formation history in <j4] 
and in Figure [5] The fraction of mass in subhalos 
and the central dominance are both highly corre- 
lated with formation time, late-time accretion rate, 
and concentration. These correlations have impor- 
tant implications for interpreting lensing-selected 
and X-ray selected clusters. 

4. A fossil cluster: Our sample i nclud es a peculiar 
outlier, Halo 572 (presented in ^4.5), with excep- 



tionally high formation redshift, concentration, and 
central dominance. It also has exceptionally low 
late-time mass accretion rate, subhalo number, and 
subhalo mass fraction. This finding indicates that 
halos of distinct formation history are likely to be 
distinguishable observationally, if stringent selection 
criteria are used. 

5. Impact of tidal stripping on the occupation 
number of subhalos: In |5j we have demon- 
strated that the subhalo number, when selected 
using Vpk (a more observationally relevant prop- 
erty), does not correlate with formation time. This 
is in contrast to the result, shown previously and 
confirmed here, that early-forming halos have fewer 
subhalos when selected with vq. We demonstrate 
that the correlation with the number of subhalos 
selected with vo can be attributed to subhalo strip- 
ping and insufficient resolution and is thus largely 
overestimated for cluster satellite galaxies. This 
finding implies that the assumption that halo occu- 
pation number is independent of formation time at 
fixed mass is likely to be a good one for luminosity 
or stellar mass selected samples, and that the for- 
mation history of clusters is unlikely to be directly 
inferred from the number of their satellite galaxies. 

In a forthcoming paper (Wu et al., in preparation), 
we will address the issue of the completeness of subhalo 
populations in detail by comparing simulations of differ- 
ent resolutions directly with observations. We will also 
investigate the impact of completeness on the measured 
velocity dispersion of subhalos. 

The lack of a correlation between halo occupation num- 
ber and formation time provides support for halo occu- 
pation models that depend only on halo mass, when the 
galaxy samples are selected by stellar mass or luminosity. 
However, our findings indicate that the halo occupation 
is likely a function of selection; the halo occupation of 
color-selected samples may depend on formation time. 
This appears to be consistent with some observational 
studies; e.g., the l umin osity-selected samples studied in 
Tinker & Conroy (2009) did not show evidence for trends 
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with formation time, while evidence for this dependence 
in samples selected by star formation rate was presented 



by Tinker et al. (|2012j). 

Although we nave found that the formation history of 
clusters does not manifest itself in the number of galaxies, 
correlation with formation time still exists for subhalo 
mass fraction, central dominance, and halo concentration. 
These correlations are potenti ally observable in tar geted 
lensing programs like CLASH ( Postman et al.|2011 ), opti- 
cal follow-up pr ograms for clusters detected by th e South 
Pole Telescope ( |High et al.|2012||Song et al|2012[ ), as well 



as the rece nt len sing mass calibration tor X-ray selected 



clusters by von der Linden et al. (2012). These properties 



together can indicate a system's state of relaxedness and 
can potentially be combined to reduce the scatter in the 
observable-mass relation for multi- wavelength surveys. 

Finally, the dependence of subhalo statistics on the 
selection method could potentially impact the prediction 
of galaxy clustering based on the halo model. For exam- 
ple, it is common to assume that the galaxy number is 
described by a Poisson distribution, or that their spatial 
distribution and velocities follow those of the dark matter 
particles (e.g. 



Zehavi et al. 2011 Cacciato et al.|2012|. As 



we have shown, these assumptions depend on the specific 
subhalo selection applied and on the simulation resolution 
and are still uncertain. Therefore, these uncertainties will 
potentially limit the accuracy with which we can predict 
the small-scale clustering and hence our ability to use it 
to infer cosmological parameters. 
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